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Abstract 

Determinations of the 3 He concentrations in Galactic matter serve to impose interesting 
and important constraints both on cosmological models and on models of Galactic chemical 
evolution. At present, observations of 3 He in the solar system and in the interstellar medium 
today suggest that the 3 He abundance has not increased significantly over the history of the 
Galaxy, while theoretical models of Galactic chemical evolution (utilizing current nucleo- 
synthesis yields from stellar evolution and supernova models) predict a rather substantial 
increase in 3 He. We consider the possibility that the solar 3 He abundance may have been 
affected by stellar processing in the solar neighborhood prior to the formation of the solar 



system. Such a discrepancy between solar abundances and average galactic abundances by 
as much as a factor of two, may be evidenced by several isotopic anomalies. Local destruc- 
tion of 3 He by a similar amount could serve to help to reconcile the expected increase in the 
3 He abundance predicted by models of galactic chemical evolution. We find however, that 
the production of heavier elements, such as oxygen, places a strong constraint on the degree 
of 3 He destruction. We also explore the implications of both alternative models of Galactic 
chemical evolution and the stellar yields for 3 He in low mass stars, which can explain the 
history of the 3 He concentration in the Galaxy. 



1 Introduction 



There is an inherent difficulty associated with the utilization of the observed abundances of 
D and 3 He to predict their primordial values. Namely, the connection between the primordial 
abundances of D and 3 He and their solar or present-day values depends sensitively on models 
of galactic chemical evolution. In principle, measurements of D in quasar absorption systems 
could dramatically help us to bridge this gap, by providing directly the primordial abundance 
of D and hence the value of the baryon-to-photon ratio, 77, from big bang nucleosynthesis 
(Walker et al. 1991). However, recent measurements of this kind (Carswell et al. 1994; 
Songaila et al. 1994; Tytler & Fan 1995) must be viewed as preliminary, as the determined 
D abundance in the two absorption systems observed are not concordant with each other. 
Until such measurements yield a single consistent value for primordial D, we must continue 
to be guided by models of galactic chemical evolution. 

It has been well established that models of galactic chemical evolution, consistent with 
the constraints imposed by element abundance determinations, are capable of destroying 
significant amounts of deuterium (Truran & Cameron 1971; Gry et al. 1984; Delbourgo- 
Salvador et al. 1985; Vangioni-Flam and Audouze 1988; Vangioni-Flam, Olive & Prant- 
zos 1994; Vangioni-Flam & Casse, 1995). However, as was recently discussed in Olive et 
al. (1995), the problem rests not with the destruction of D, but rather with the production 
of 3 He. Though 3 He is partially destroyed in massive stars, 3 He production in low mass 
stars generally leads to a net increase in the 3 He abundance over the evolutionary history 
of the galaxy. Observations of large 3 He enhancements in planetary nebulae (Rood, Bania 
& Wilson, 1992; Rood, Bania, Wilson & Balser, 1995) support the conclusion that 3 He is 
indeed produced in low mass stars. An excess of 3 He is difficult to avoid if low mass stars 
are strong producers of this isotope, as indicated by the calculations of Iben and Truran 
(1978) and more recently Vassiliadis and Wood (1993) and Weiss et al. (1995). However, 
recent models with increased mixing have been calculated, which both bring the carbon and 
oxygen isotopic ratios to their observed level in red giants and lead to a net destruction of 
3 He (Charbonnel 1994; Wasserburg et al. 1995; Hogan 1995). At present, these models are 
quite preliminary, and it is premature to draw a firm conclusion. However, if their results 
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are confirmed, 3 He will be significantly less problematic, unless the primordial D abundance 
is as high as observed by Songaila et al. (1994) and Carswell et al. (1994); in this case we 
would still require 3 He destruction factors in excess of what current calculations bear out 
(Olive et al. 1995). 

Abundance determinations of 3 He at the time of the formation of the solar system seem 
to indicate that the solar 3 He abundance is very close to that of the primordial abundance. 
Here, we will examine in some detail the possibility that the 3 He abundance in the solar 
system may be depleted with respect to galactic averages at the time of its formation. We 
find, however, that the abundances of the heavier elements, most notably oxygen and neon, 
impose a strong constraint on the degree of depletion of 3 He. We are therefore left with the 
following possibilities: either the initial deuterium abundance is low, D/H ~ 3 x 10 -5 (we 
will discuss rough lower limits imposed by models of chemical evolution); or more dramatic 
changes are required in models of chemical evolution, which have the effect of maintaining a 
rather flat evolution of 3 He with time (we will show an example of this type of model below); 
or the stellar yields of 3 He in low mass stars are lower than previously thought. 

Because deuterium is converted to 3 He in the pre-main-sequence phase of stellar evolution, 
models without a significantly depressed initial D abundance are subject to problems with 
3 He production. If the initial abundance of D is rather low (D < 3x 10~ 5 ), the present day 
3 He abundance found in standard models of galactic chemical evolution are not excessive 
and are consistent with the observed range of 1 — 5 x 10~ 5 (Balser et al. 1994). However, 
even in these cases, there appears to be a problem with the abundance of 3 He as measured 
in meteorites, giving the pre-solar value of 3 He. That is, on the basis of chemical evolution 
models, we expect more 3 He than is observed in the solar system (Olive et al. 1995; Galli et 
al. 1995; Tosi et al. 1995). In this paper, we will thus consider in turn, the possibilities that 
3 He in the solar system has been depleted with respect to the galactic average; the accuracy 
of the measurement of solar 3 He; and what we can learn from galactic chemical evolution. 
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2 Solar Depletion of 3 He? 



In principle, it is possible that the abundance of 3 He at the time of the formation of the 
solar system does not reflect the galactic average at that time. It is necessary to consider the 
degree to which element abundances in the solar system were affected by the explosions of 
supernovae in the solar neighborhood, immediately prior to the formation of the solar system 
(Reeves 1978; Olive & Schramm 1982). This may be evidenced by several "anomalous" 
isotopes (of Carbon, Oxygen and Neon) seen in cosmic rays. However, late contributions to 
the abundances of 16 O and 20 Ne, which are produced solely in Type II supernovae, may render 
the solar system isotopic ratios of these elements anomalous. That is, the solar abundances 
may not represent the true average galactic abundance. For the oxygen and neon isotopes, 
these differences may be as large as a factor of two. We now examine the possibility that 
the solar 3 He abundance may also not be representative of the true galactic abundance. 

From the observed anomalies in 26 Al and 107 Pd (Lee 1979; and references therein) and 
more recently 41 Ca, and the short time scales associated with their half-lives (~ 10 6 years), 
the element abundances in our solar system were probably affected by at least one supernova 
within that time period prior to formation. Even a single supernova explosion in a star 
forming region can have dramatic consequences on the element abundances of that region. 
As was argued by Olive & Schramm (1982), a handful of the first few supernovae in an 
early OB association, is capable of producing nearly the entire observed solar abundance of 
16 and 20 Ne. Thus we would expect that the solar isotopic ratios such as 17 ' 18 0/ 16 and 
22 Ne/ 20 Ne may be diluted with respect to the galactic average. It is therefore of interest to 
question whether or not the abundance of 3 He (which would be depleted in the ejecta of the 
first few supernovae in an association) is comparably depleted. 

To deplete 3 He, we must require that a significant amount of material in the solar neigh- 
borhood underwent stellar processing prior to the formation of the solar system. Let us 
suppose that a fraction / of the total initial gas of the association went into stars prior to 
the solar epoch. It is then reasonable to assume that a fraction ~ 0.1/ of the gas went into 
stars with masses greater than 10 M Q . (For example, with a Salpeter (1955) initial mass 
function (IMF) <f>(m) oc m~ 2 - 35 between, 0.1 and 100 M , the fraction is 12%; for a Scalo 
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(1986) mass function it may range from 5-15% depending on the star formation rate (SFR)). 
If we denote by X* the mass fraction of some heavy element (such as O) ejected by massive 
stars, the total mass fraction of the element after the explosions of stars more massive than 
10 M & which thus determines (and must be less than) the solar abundance, is then given by 

0.1/X. + (!-/)*< 

Xf ~ (i - /) + 0.1/ <X& (1) 

where Xj is the initial mass fraction of the element and we have assumed that 0.9/ of the 
initial gas mass is still locked in stars. To maximize our estimate of 3 He destruction (as this 
will maximize our estimate for /), we can assume that X, = 0. Solving for /, we have, 

f < ^® ( 2 ) 

J O.1X* + O.9X 1 ; 

For oxygen, X* ~ 0.1 and X Q ~ 0.01, so that / <, 1/2. Thus, we can cycle no more than 
1/2 the mass of the association through stars prior to the formation of the solar system. 

Although a significant amount of gas may be cycled through stars, only a small fraction 
of 3 He depleted gas can be released back into the association. If we take X to represent 3 He 
in Eq. (P, and now take X* = (which assumes that 3 He is totally destroyed in massive 
stars), and Xj to be the primordial 3 He mass fraction, we find 

x > = T^f ~ sx < < 3 > 

This indicates that only about 10% of the initial 3 He can be destroyed, even though changes 
in the heavy element abundances occur at a level of a factor of 2. 

It is possible, of course, to further deplete 3 He in the gas which forms the solar system 
at the expense of excessive metal production. Such overproduction of metals can perhaps 
be reconciled with 3 He depletion, if the heavy elements could somehow be expelled from 
the solar neighborhood. As Lattimer, Schramm & Grossman (1977) pointed out, the bulk 
of the heavy element ejecta from supernovae can rapidly form into dust grains. These dust 
grains can behave like explosive "shrapnel" and penetrate regions exterior to the association. 
This might allow the association itself and hence the solar system, to fail to show a large 
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heavy element excess, even though the total heavy element enrichment would be part of the 
integrated galactic enrichment. This assumes, of course, that the entire association region is 
not totally disrupted by the supernovae explosion prior to the formation of the solar system. 
However, as was shown in Olive & Schramm (1982), a significant amount of oxygen and 
neon is produced which should not be trapped in grains. Because of the behaviors of these 
elements, we believe that is unlikely that more than about 10% of the 3 He present in the 
association could be destroyed before the formation of the solar system. 

An obvious recourse to resolving the problem of the overproduction of 3 He at the solar 
epoch, is to question the measurement of the solar 3 He abundance. In the next section, we 
examine the observational data on the solar abundances of D and 3 He. 

3 D and 3 He in Pre-solar Nebulae 

Because the crucial data in attempts to estimate the primordial D/H and 3 He/H values 
come from solar system measurements, it is useful to examine critically the origin of these 
abundances and to attempt to provide an accurate estimate of their uncertainties. 

The determination of the solar abundances of D and 3 He involves 3 He/ 4 He measurements 
both in meteorites and also directly in the solar wind. Direct D/H measurements are irrele- 
vant for the Sun, since D is completely burned to 3 He in the solar convective zone. Moreover, 
D/H measurements are difficult to interpret in planetary bodies (Earth, Jupiter, etc.); be- 
cause D preferentially enters molecules relative to H, abundance determinations thus require 
a knowledge of complex chemical fractionation histories. However, because essentially all pri- 
mordial D has been burned to 3 He in the solar convective zone and because the convective 
zone is not hot enough to burn 3 He, the solar wind measurement of 3 He/ 4 He provides a mea- 
surement of the pre-solar abundance of D t H ^ e I o by number. Solar wind measurements, made 
using foil collectors during Apollo lunar missions, yielded values for 3 He/ 4 He ranging from 4 
to 5.5 xlO" 4 . Geiss and Reeves (1972) and Bochsler & Geiss (1989) (see also Geiss (1993) for 
a recent review) argue that the variation can be corrected for, and that the best solar wind 
ratio is 3 He/ 4 He| sw = 4.1 ± 1 x 1CT 4 (where the error is statistical). This is in good agree- 
ment with the low temperature component emitted by carbonaceous chondrites in step-wise 
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heating experiments (Black 1972; Weiler et al. 1991), for which 3 He/ 4 He| sw ~ 4.5 ± 1 x 10~ 4 , 
and also with the ISEE-3 solar wind data (Coplan et al. , 1984), which yields 4.4 x 10~ 4 . 
However, some fractionation in all the solar wind 3 He/ 4 He measurements cannot be ex- 
cluded, which would add an additional systematic error to the above value. The extreme 
value for 3 He/ 4 He observed for the Apollo solar wind measurement of 5.5 x 10~ 4 can not be 
excluded as a central value, hence a systematic uncertainty of 1.4 x 10~ 4 for 3 He/ 4 He cannot 
be excluded at present. The most recent measurement of 3 He/ 4 He in the solar wind from 
the over-the-solar-pole measurements made with the SWICS instrument on the ULYSSES 
spacecraft gives 3 He/ 4 He = (4.4 ± 0.4) x 1(T 4 (Bodmer et al. 1995). 

The pre-solar 3 He/ 4 He| ratio is thought to be best measured in meteorites. Initially, 
Black (1971) proposed that the high temperature component emitted by step-wise heating 
experiments using carbonaceous chondrites (see also Eberhardt 1974) was the primordial 
component. 3 He/ 4 He ~ 1.5 x 10 -4 . However, Weiler et al. (1991) have argued that this high 
temperature component is dominated by gas trapped in pre-solar grains (diamonds) which 
formed in locations far removed from the solar system. Weiler et al. propose that another gas 
component known as "Q" is a better candidate for the primordial component. Fortunately, 
the difference in 3 He/ 4 He between the high T carbonaceous chondrite component and Q is 
relatively small 



However, a potential interpretational (systematic) error persists, since neither Q nor dia- 
monds nor the high T carbonaceous chondrite component has been unequivocally proven to 
represent 3 He/ 4 He| . Taking 3 He/ 4 He| g as 3 He/ 4 He| , but allowing for systematics to include 
the range of relevant meteoritic 3 He/ 4 He values, yields 3 He/ 4 He| = 1.6 ± 0.04 ± 0.3 x 10~ 4 . 
The pre-solar D is estimated by subtracting 3 He/ 4 He| from the SWICS solar wind value. To 
convert to ratios relative to hydrogen requires multiplying by the number ratio of 4 He/H| , 
which is estimated to be 0.09 ± 0.01 (note, this is 10% lower than that used by Geiss 1993) 
from the best fit solar model Y = 0.27 (Turck-Chieze et al. 1988; Bahcall and Pinnsoneault 
1992) with metallicity Z = 0.02. This yields 



3 He/ 4 He\ Q = 1.6 ± 0.04 x 10 



-4 



(4) 



D+ 3 He 
H 



4.1 ± 0.6 ± 1.4 x 10 



,-5 



(5) 
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and 



3 He 



1.5 ±0.2 ±0.3 x 10 



,-5 



(6) 



H 



and thus 



D 



,-5 



(7) 



2.6 ± 0.6 ± 1.4 x 10 



He 



This latter number is in reasonable agreement with the ^ = 1 — 3 x 10~ 5 ratio measured 
in Jupiter (Smith, Scherpp & Barnes, 1989). Although planetary D ratios are subject to 
chemical fractionation, this is minimized for HD on Jupiter, since the bulk of the deuterium 
and hydrogen is in HD and H2 there. However, molecular line blanketing does still allow for 
significant systematic errors. For this reason, Jupiter is still not the best source for a solar 
system D determination, but it does provide a consistency check. 

4 Chemical Evolution 

The solar system abundance of 3 He is thus seen to be approximately a factor of two lower 
than predicted by even the more optimistic models of Galactic chemical evolution, which 
tend to yield abundance ratios at least as high as 3 x 10~ 5 for 3 He/H, when 3 He production in 
lower mass stars is included (Olive et al. 1995). In what follows, we will look at three different 
approaches to resolving the problem of excess solar 3 He. We first consider possibilities for 
which the primordial value of D/H is low. A low initial D/H lowers 3 He/H, as there is less 
D to be converted to 3 He in the pre-main-sequence evolution of stars. However, as we will 
show, one can not take arbitrarily low values of D/H (of course D/H is always bounded from 
below by the ISM measurements of D/H yielding D/H = 1.6 ± 0.09tg;f (Linsky et al. 1993, 
1995)), since some amount of deuterium destruction necessarily accompanies the production 
of heavy elements in the galaxy. We then consider "higher" values of D/H, which require some 
dramatic changes to simple models of chemical evolution, such as an increased production 
of massive stars in the early galaxy as well as metal enriched outflow. We will also examine 
some remaining alternatives regarding the stellar production of 3 He. Note however, that 
there may be a quite disturbing dispersion of D/H in the local ISM which would complicate 
the analysis (Vidal-Madjar 1991; Ferlet 1992, Linsky private communication) 
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As was noted earlier, the questions concerning high verses low D/H may become moot, if 
the determinations of primordial D/H in quasar absorption systems yield a single consistent 
value. To date there are three measurements of D/H in quasar absorption systems. Two 
(in the same system) yield a high value for D/H ps 1.9 — 2.5 x 10~ 4 (Carswell et al. 1994; 
Songaila et al. 1994), while the third (in a different system yields a significantly lower value, 
D/H 1 — 2 x 1CT 5 . It is clear that, on the basis of these measurements, we can not with 
confidence claim any knowledge of the primordial abundance of deuterium. Indeed, it has 
been argued ( Levshakov & Takahara 1995) that measurements of this type may not be able 
to determine D/H to better than an order of magnitude. In other words, they would expect 
a large dispersion in the observational data. Is this what we are seeing? 

Interestingly enough, the two values for D/H identified above are in some respects both 
beneficial and detrimental to big bang nucleosynthesis. The high value of D/H corresponds to 
a value for the baryon-to-photon ratio 77 ~ 1.5 x 10~ 10 (Walker et al. , 1991). Consequences of 
this high D/H were recently discussed in Vangioni-Flam & Casse (1995). With regard to the 
other light elements produced in big bang nucleosynthesis, the low value for 77 corresponds 
to a 4 He mass fraction Yp ~ 0.23, which is in remarkable agreement with what one expects 
from the data on 4 He from extragalactic H11 regions (Olive & Steigman 1995; Olive & Scully 
1995). 7 Li/H is predicted to be around 2 x 10~ 10 which is also compatible within errors, with 
recent data (Molaro et al. 1995). The problem occurs with the evolution of 3 He, when 3 He 
production is included (note that models of chemical evolution can be constructed which can 
account for the necessary D/H destruction in this case). In Olive et al. (1995), it was found 
that the abundance of 3 He at the time of solar system formation could be high by as much 
as a factor 10. Even in the absence of 3 He production, it was found that massive stars were 
required to destroy at least 90% of their initial D + 3 He, in order to reproduce the solar and 
ISM values of 3 He. This amount of destruction is excessive, even for the most massive stars 
(Dearborn, Schramm & Steigman, 1986). 

On the other hand, the low value of D/H between 1 and 2 xl0~ 5 corresponds to a 
value of 77 ~ 7 — 9 x 10~ 10 . In contrast to the high D/H case, we would expect a much 
milder problem with 3 He (to be discussed below). However, now the 4 He mass fraction is 
predicted to be Y P > 0.249, a value larger than most of the 4 He measurements (Pagel et 
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al. , 1992; Skillman et al. 1995) in extragalactic Hn regions, which already contain some 
non-primordial 4 He. (However, again, possible systematic errors can not be excluded Copi 
et al. , 1995a; Sasselov & Goldwirth, 1995.) In addition, 7 Li/H is expected to be > 5 x 1CT 10 
requiring a significant amount of 7 Li depletion, contrary to what one expects (Steigman et 
al. 1993) from the positive measurements of 6 Li in halo stars (Smith Lambert & Nissen, 
1992; Hobbs & Thorburn, 1994). Furthermore, as we will next show, a minimal amount of D 
destruction is demanded for consistency with the observed level of heavy element production 
in the Galaxy. A completely flat evolution for D is probably excluded on these grounds. 

The classical constraints on galactic evolution are characterized by varying degrees of 
stringency. Among these, the trends in [Fe/H] with time are easily satisfied, since the age- 
metallicity relation suffers from a large dispersion over the observed age range (Edvardsson 
et al. 1993; Nissen 1995). The [O/Fe] vs [Fe/H] relationship is mainly sensitive to the 
stellar yields and not to the different histories of star formation (assuming a constant IMF). 
The metallicity distribution of disk stars is far from being definitely established. Indeed, 
much work is needed before a clear picture of the metallicity distribution can be reached 
(e.g. Olsen 1994; Cayrel, private communication). Information on metallicities, ages and 
kinematics, with the same high accuracy as obtained by Edvarsson et al. (1993), is needed 
for a much larger stellar sample. Moreover, Grenon (1989, 1990) remarks that the radial 
migration of stars in the Galaxy can blur the local metallicity distribution. 

Other global characteristics which should be considered are the gas fraction, a, the over- 
all metallicity, Z, and individual abundance ratios (Fe/H, O/H,...) at solar birth and in 
the present ISM. To the list of constraints, we must also add the D/H and 3 He/H ratios 
at solar birth and at present time, in relation to the primordial value. Indeed, since pri- 
mordial nucleosynthesis is much more constrained than galactic evolution, it is reasonable 
to harmonize the second to the first, and not the contrary (as has sometimes been done 
recently). 

Many models have been proposed to follow the chemical evolution of the Milky Way, 
invoking, for example, a prompt initial enrichment ( Truran & Cameron 1971), infall of 
primordial material (Timmes et al. 1995; Fields 1995), metal enriched infall originating from 
the halo (Ostriker and Thuan 1975), and early massive star formation ((Larson 1986; Wyse 
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and Silk 1987). Studies of galactic chemical evolution remain in their infancy, however, since 
we do not yet have good theories of galaxy formation and star formation. It would be unwise, 
for the sake of simplicity, to limit the investigation to "classical" models under the pretext 
that they have been widely used. In effect, if the high primordial D/H ratio is confirmed, 
special models leading to a strong D destruction avoiding overproduction of 3 He and Z will 
be required. 

An alternative way of looking at variations from the galactic mean has been carried 
out by Copi, Schramm & Turner (1995b) looking at the stochastic variations from galactic 
evolution models. Their conclusions concerning the allowed range of primordial D and 3 He 
are similar to, and compatible with, those we discuss here. 



We will first explore the possible consequences of a very low primordial value of D/H and 
examine the extent to which a low D/H could explain the apparent flatness of the 3 He/H 
evolution in the Galaxy. We begin by estimating the minimum possible amount of D/H 
destruction. In simplified models of galactic chemical evolution, it is possible to derive some 
analytic relations between abundances, yields, the gas fraction, and the IMF, if one assumes 
the instantaneous recycling approximation (that is, that the enriched mass that is ultimately 
to be ejected from a star is incorporated into the ISM at the time of formation of the star, 
in contrast to its appropriate delayed entry at the end of the star's lifetime). Indeed, the 
degree to which deuterium is destroyed can be expressed simply by (Ostriker and Tinsley 



In (^|), Mi is the main-sequence turnoff mass (normally a function of time), M sup is the upper 
mass limit for star formation, and M rem is the remnant mass. 



4.1 Low D/H 
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(9) 
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It is also possible to express the metallicity in terms of the gas mass fraction and the 
yields of metals in stars (Searle & Sargent 1972) 

Pz 



(l-R) 

where 



In a" 1 (10) 



Msup My 

-jM0M)dM (11) 
I Mi M 

and Mz/M is the mass fraction ejected in metals. Equations (||) and ( |10"1) can be combined 
yielding 



D 



(12) 



As one can see from Eq. (|8]), a low primordial value for D/H, will require a small return 
fraction, R. In principle, one can easily adjust the IMF to yield a small value for R. However, 
because of the similarity in the definitions of R and Pz-, their ratio is almost independent of 
the details of the IMF. Thus D/D p near unity, implies a metallicity much less than solar. 

The interdependence between deuterium and metallicity can be seen in Figure 1. In 
order to reach solar metallicity at the time the solar system formed, we require a deuterium 
destruction factor of at least 1.6, implying that D/H p >, 2.5 x 10~ 5 . We note that this factor 
is somewhat dependent upon the assumed yields for the heavier elements. For example, this 
limit was obtained using the stellar yields of Woosley & Weaver (1993), whereas had we used 
the yields of Maeder (1992), which allow for more heavy element production in the mass range 
from 9-11 M , the minimum destruction factor could be lowered to about 1.3. It is worth 
noting that, beyond the uncertainties of the yields which are essentially related to those 
associated with the 12 C(a,7) 16 reaction rate, the lower mass limit of the stellar progenitor 
of the Type II supernovae which synthesize the heavy elements is influential because of the 
preference in the IMF towards lower mass stars. Indeed the limit is greatly increased as 
the slope of the IMF is decreased. All of these effects can be seen in Figure 1, where we 
have plotted (for various choices of the parameters which govern the SFR) the metallicity at 
the solar epoch in units of solar metallicity, Z/Z , as a function of the ratio of the present 
deuterium abundance to the primordial one, thus indicating the total deuterium destruction 
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factor for a variety of galactic evolution models. We have chosen a SFR proportional to the 
mass in gas, and an IMF, 4>(m) oc m -2 ' 7 ; shown here by the upper set of points denoted 
by circles (yields from Maeder (1992)) and crosses (yields from Woosley & Weaver (1993)). 
For the lower set of points, a steeper IMF, (f>(m) oc m -3 was chosen. In each case, the lower 
mass limit of the IMF was 0.4 M (lowering this choice to 0.1 M would further lower the 
curves). It is important to note that, when 3 He production is taken into account, even the 
modest deuterium destruction factor (of 1.6), yields an overproduction by about a factor of 
2 in solar 3 He. 

4.2 Higher D/H 

In this section, we will consider an alternative to low primordial D/H and rely on more 
distinctive models of galactic chemical evolution to resolve the problem concerning the solar 
3 He abundance. As we have stated earlier, the choice of a higher value for primordial 
D/H alleviates some of the pressure in matching the BBN calculations to the observational 
determinations of 4 He and 7 Li. Clearly, the higher the value we choose for primordial D/H, 
the more difficult it will be to keep 3 He under control. We choose specifically the value 
D/H| p = 7.5 x 10" 5 which corresponds roughly to the 7 Li trough and is in modest agreement 
with 4 He (at the 2 a plus systematics level). 

The models we consider below specifically involve mass outflow. Open galactic models 
have been considered in the past (Tinsley 1980, Tosi 1988), but infall has been invoked 
more often than outflow. Formally, the two reverse processes are included in the general 
formalism of chemical galactic evolution (Tinsley 1980) and cosmochronology (Cowan et 
al. 1989). There is clear evidence for galactic winds in external galaxies, even for spirals 
(Wang et al. 1995), and particularly those experiencing bursts of star formation, whereas 
evidence for significant infall of extragalactic matter are meager (Murphy et al. 1995). Of 
course this reflects only constraints arising from the present state of the solar vicinity, and 
proves nothing about the early galaxy. 

Outflow has its own merit: as we will see below, it will help to explain very high destruc- 
tion factors of D (Vangioni-Flam & Casse 1995) while, if necessary, avoiding, an overproduc- 
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tion of metals. At the same time it could reduce the rise of the 3 He/H ratio. De Young & 
Heckman (1994) proposed that energy from supernova explosions and stellar driven winds 
results in blowing portions of the ISM containing enriched material out of the galaxy. Differ- 
ent kinds of outflows can be imagined, which vary considerably in their durations, intensities 
and compositions. For our present purposes, it is sufficient to distinguish whether the out- 
flowing matter consists solely of the ejecta of massive stars, or rather whether it is composed 
of normal ISM material being blown out by supernova driven winds. In Casse et al. (1995), 
we will return to these distinctions in greater detail. 

The three specific ingredients that must be added to canonical galactic evolutionary 
models in order to obtain significant D destruction without the overproduction of 3 He are: 
(i) an early phase of massive star formation, which presents the advantage of destroying D 
and 3 He rapidly; (ii) a galactic wind related to the corresponding SNII rate, which limits 
the rise of Z and 3 He, leading to an even more pronounced decrease of D; and (iii) possible 
modifications of stellar models, leading to an efficient destruction of 3 He, especially in low 
mass stars. 

One way we have found in which the solar value of 3 He may be lowered is to assume that 
the IMF prior to the formation of the solar system was skewed more toward massive star 
formation. The presence of fewer lower mass stars reduces 3 He production, while more mas- 
sive stars ultimately return only a fraction of the 3 He present during the pre-main sequence 
phase to the ISM. We therefore consider models which begin with an IMF favoring more 
massive stars early on galactic history but resemble a more normal IMF at later times. 

The problem that we immediately encounter is that the emphasis on more massive stars 
results in an overproduction of heavy elements, such as 16 0. We have found that the (closed 
box) models which are most successfully in keeping 3 He flat, while destroying enough D, also 
overproduce 16 by a factor of ~10. This problem is alleviated by including outflow. Indeed, 
McCray & Snow (1979) have shown that supernovae can generate "chimneys," which can 
directly transport much of the heavy element rich supernova debris out of the galaxy. We 
have therefore included "enriched" (relative to the ISM) outflow in our models, both to help 
solve the heavy element overproduction problem and to obtain a flatter 3 He/H evolution. 
In order to simulate this effect, we have incorporated outflow into our models at a rate 
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proportional to the rate of ejection of materials from supernovae. 

Since massive stars can lose large amounts of 3 He depleted outer material via winds 
before they explode, it is certainly possible for them to deplete 3 He in their surrounding ISM 
material and eject their metals out of the Galaxy. We allow the outer (hydrogen) envelope of 
the star which is deficient in 3 He to return via winds to the ISM. Then, in order to maximize 
the possible effect of an outflow which is tied to the ejecta of exploding stars (M > 8-10 
M Q ), a fraction of the core is then expelled from the Galaxy. Such models provide a natural 
way to understand the heavy element abundances in the X-ray gas observed in clusters of 
galaxies (typically, the heavy element enriched outflow is produced by elliptic galaxies, (see 
e.g. Elbaz, Arnaud, Vangioni-Flam 1995)). It might also be noted that early expulsion of 
metal rich supernova ejecta is even easier in merger models, where the early galactic building 
blocks have a lower mass. Due to the epoch of more massive star formation at earlier times, 
D is generally very efficiently destroyed in these models. 

We should note at this point that our assumptions concerning winds from massive stars 
prior to the supernova stage may be inappropriate at early galactic epochs. The rate of mass 
loss is generally expected to be dependent upon the initial metallicity of the star (Maeder, 
Lequeux, & Azzopardi 1980; Maeder & Meynet 1994). Expectations from theoretical studies 
are generally consistent with e.g. trends in the frequency of Wolf-Rayet stars, as inferred from 
studies of the Magellanic Clouds (Massey et al. 1995). This suggests that the fraction of the 
3 He depleted outer envelopes of massive stars that is returned to the ISM via winds (prior 
to supernova-triggered mass ejection) in low metallicity populations may be significantly 
reduced. We stress however, that our aim here is to see how efficiently the evolution of the 
3 He abundance can be held relatively flat over the history of the Galaxy. As we will see, we 
find only modest success despite rather poignant assumptions. 

An obvious observational constraint on our choice of an IMF that is skewed toward more 
massive star production at early times, (j>(m) oc m _ ( L25+ °/° ), is provided by its consistency 
with the present day IMF that results from our model. Fig. 2 shows a comparison of the 
observed and modeled present day IMF. The observed values are taken from Scalo (1986) 
and are in good agreement with our model for the more massive stars. This is as expected, 
since the more numerous massive stars formed early on have long since died out. 
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It is be useful to compare the results we have obtained here with those of our previous 
study, in which we considered a more standard model. In model 1 of Olive et al. (1995), 
we chose a SFR, ip = 0.25a, with an IMF, (f>(m) oc m~ 2 ' 7 between 0.4 and 100 M . The 
3 He abundance at the solar formation epoch (taken to be at t = 9.4 Gyr) was 3 He/H = 
5.7 xlO -5 , rising to 8.8 xl0~ 5 today. Infall was not included. For a primordial ratio D/H 
= 7.5 x 10~ 5 , the results for D/H and 3 He/H as a function of time, for a model with an 
IMF which favors massive stars early, and contains enriched outflow in which the rate is 
proportional to the ejection rate, is compared with model 1 of Olive et al. (1995) in Figure 
3. The outflowing gas contains only material below the outer envelope, while the latter 3 He 
depleted material is returned to the ISM. In this model, the SFR is ip = 0.26M gas , and the 
fraction of outflowing gas is 90% of the supernova ejecta. The IMF is now extended down 
to 0.1 M , to help keep the metallicity and gas mass fraction reasonably low. We view this 
as a rather extreme model, in which a considerable amount of enriched material has been 
expelled from the galaxy. Indeed, we impose a limit, arising from the observed metallicity 
of hot X-ray gas in clusters, on the amount of metals expelled by outflow to be less than 
20 times the amount of metals in the galaxy. This imposes a constraint on the fraction of 
outflowing gas (90% in this case). 

As one can see from Fig. 3, our present model, which is based on an IMF skewed 
towards massive stars early on and contains enriched outflow, reduces the abundance of 3 He 
by a factor of about 2, relative to the standard case with a normal IMF and no outflow. 
Parameters of the model have been chosen such that the degree of deuterium destruction is 
comparable (and agrees with the data) in the two cases. However, although the present 3 He 
abundance is acceptable 3 He/H| Q = 5.1 x 10 -5 , the solar abundance is still high by a factor 
of slightly over two; 3 He/H| = 3.7 x 10~ 5 . While this represents a definite improvement, 
it can not be regarded as a solution to the problem. Although it appears from Figure 3, 
better agreement with the solar data is possible if one assumes a lower time for the formation 
of the solar system, the model must be adjusted to destroy D on a faster time scale. For 
example, with if) = .34M gas the evolution of deuterium matches the solar (and present-day) 
observations at t — 6 Gyr (corresponding to an age of the Galaxy of 10.6 Gyr), but now 
3 He/H| ~ 3.1 x 10~ 5 and the present abundance is 4.9 xlO -5 ; a further improvement, but 
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solar 3 He is still too high. 

Of course as is well known, the problem concerning 3 He is also alleviated somewhat by 
going to higher values of rj. In model 3 of Olive et al. (1995), we assumed a primordial 
abundance of deuterium of D/H = 3.5 xlO -5 . In this case 3 He/H| = 3.4 x 10~ 5 , an 
overproduction by a factor greater than two. The present 3 He was also slightly high, 3 He/H 
= 6 x 10~ 5 . In models with outflow as described above, these numbers are reduced to 
3 He/H| = 2.3 x 1(T 5 and 3 He/H = 3.2 x 1(T 5 today. 

Before, we move on, we wish to stress that the problems concerning 3 He that we are 
discussing here, prevail only because we are including the production of 3 He in low mass 
stars. When such production is ignored there is no problem in matching the solar and ISM 
data for both D and 3 He in models of these types as was shown by Vangioni-Flam, Olive 
& Prantzos (1994). The crises in big bang nucleosynthesis claimed by Hata et al. (1995) is 
only a crises because of the limit on the degree of 3 He destruction they allowed. Although 
the final 3 He abundance in a given star relative to the initial D + 3 He abundance, usually 
called (?3, is always larger than 0.25 as assumed by Hata et al. , even simple models such as 
the type considered here (without outflow) and in Vangioni-Flam, Olive & Prantzos (1994) 
have an effective g% which is lower than 0.25 vitiating the purported crises. 

4.3 Alternatives 

A critical consideration with regard to the establishment of any realistic constraints on 
cosmological D and D + 3 He is that associated with 3 He production in low mass stars. 
Essentially all early estimates of D and 3 He constraints on cosmology (see, e.g., Truran & 
Cameron 1971; Rood, Steigman, & Tinsley 1976) were based upon the stellar evolution 
models of Iben (1967 a,b), for which analytical fits to the detailed model characteristics were 
subsequently provided by Iben & Truran (1978). The problem of 3 He then is simply the fact 
that, with the use of the Iben & Truran (1978) prescriptions, 3 He production in stars in the 
mass range ~ 1-3 M is sufficient to overproduce 3 He in Galactic chemical evolution models 
(Olive et al. 1995; Galli et al. 1995; Timmes & Truran 1995), relative both to the solar 
system value of 3 He and to the 3 He concentration in the ISM at the present time (Balser et 
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al. 1994). Further strong confirmation of this behavior has been provided by recent stellar 
evolution calculations (Vassiliadis & Wood 1994; Weiss et al. 1995). It would seem to be 
necessary either to utilize rather extreme assumptions regarding the history of our Galaxy 
or to identify some significant problem in stellar evolution theory. 

An interesting recent paper by Wasserburg, Boothroyd, & Sackmann (1995) has called 
attention to the fact that the long-standing problems associated with understanding both 
low 12 C/ 13 C ratios in low mass red giant branch stars and low 18 0/ 16 ratios in asymptotic 
giant branch stars can be resolved, with the assumption of the occurrence of deep circulation 
currents extending below the bottom of the standard convective envelope. A concomitant 
of this process of "cool bottom burning" is the destruction of 3 He. In particular, for the 
case of a 1 M star, their models predict that after a I s * dredge-up 3 He enhancement of a 
factor ~ 6, cool bottom processing acts to reduce the 3 He concentration by a factor ~ 10, 
yielding a net depletion of 3 He by a factor ~ 2. If this model is indeed correct, this would 
aid substantially in the problem of 3 He overproduction, with which we are so concerned in 
this paper. 

To test the effect of such a reduction in the 3 He yields, we incorporated the results of 
Wasserburg et al. (1995) by lowering the Iben & Truran (1978) yields of 3 He at 1 M Q by 
a factor of 10. For an initial deuterium abundance of 7.5 xlO -5 , this corresponds to a 
gs = 0.27. We reduced the degree to which the Iben & Truran yields were modified at higher 
masses such that, at M > 3 M , we once again were using the Iben & Truran yields. The 
results of such a reduction in model 1 of Olive et al. (1995) are shown in Figure 4. Here 
the evolution of D/H and 3 He/H are shown in model 1 with both the Iben & Truran yields 
and the reduced yields. Even in this case, there remains a mild overproduction of 3 He by 
a factor of about 2. That is, at t = 9.4 Gyr, 3 He/H = 3 xl0~ 5 . In Figure 5, we show the 
effect of the reduced 3 He yields in the model with outflow discussed in the previous section. 
Here finally, we find a value for 3 He/H at the solar epoch which is perhaps acceptable. At 
t = 9.4 Gyr, 3 He/H = 2.4 xlO -5 . A further improvement is possible by considering models 
which evolve on shorter time scales as discussed above. 

An obvious problem with the reduction in 3 He yields at low stellar masses is the ob- 
servation of high 3 He concentrations in planetary nebula ejecta (Rood et al. 1992, 1995), 
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which would seem to confirm the predictions of the more standard models for the evolution 
of low mass stars along the giant branch. It is clear that this issue must be resolved before 
a more definitive statement can be made with respect to Galactic evolution constraints on 
the primordial abundances of D and 3 He. 

A further question of interest is that concerning the composition of the matter involved in 
"mass infall," during the later stages of evolution of our Galaxy. In this context, we note that 
while we have not considered such infall models in this paper (see, e.g., Olive et al. 1995), 
they may provide plausible alternative solutions to the 3 He problem. The implications of 
infall of matter of primordial composition of the light elements D. 3 He, and 7 Li are certainly 
quite different from those of processed matter, which may generally be expected to be metal 
enriched and deuterium depleted. Infall of primordial material is generally beneficial, with 
the adoption of the lower primordial D abundance, while infall of D- and 3 He-depleted matter 
improves the situation for the case of a higher primordial D abundance. The fact that the 
nature and origin of such infalling material is presently uncertain, makes it necessary to 
treat its composition as an additional parameter. This problem is further complicated by 
the fact that it is even possible for the in falling gas both to be metal enriched and to have an 
essentially primordial composition of D and 3 He. Such could occur if, for example, the ejecta 
of the first generation of massive stars in the halo of our Galaxy were lost to the surrounding 
intergalactic medium. The ejecta of stars of M > 1OM collectively represents ~ 10 % of 
the initial mass formed into stars (e.g., for a Salpeter IMF over the range 0.1-100 M©) and 
is characterized by a metal abundance > 10 Z . Assuming ~ 10 10 M Q of early halo star 
formation would yield ~ 10 9 M of metal enriched gas ejected, which could contaminate ~ 
10 10 M to solar metallicity and yet have deuterium at a level of only ~ 0.9 its primordial 
value. 

5 Conclusions 

What can we conclude from this analysis? We have shown that, unlike the the abundances 
of some of the heavier elements such as oxygen and neon which can differ by as much as a 
factor of two locally relative to their average galactic abundance by prior supernova in the 
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solar neighborhood, the local 3 He abundance could only have been affected by at most 10%. 
It also appears that the 3 He data from a variety of sources is consistent and yields a value 
3 He/H| = 1.5 x 10 -5 for the presolar 3 He abundance. Standard models of galactic chemical 
evolution yield an excess of 3 He at the solar epoch by a factor which ranges from 2 to 12 
depending on the assumed primordial value for D/H. For models with D/H = 7.5 x 10~ 5 
initially, the factor nearly 4 excess in 3 He can be brought down to an excess of about 2 
in models which favor massive stars early on, and include the possibility for a substantial 
amount of metal enriched outflow. In such models, the solar 3 He abundance is brought down 
to nearly acceptable levels when primordial D/H < 3.5 x 10~ 5 . Finally, we considered the 
possibility that part of the problem may lie in the stellar yields of 3 He. Though it appears 
that the cut in yields suggested by Wasserburg et al. (1995) may not in itself be sufficient to 
lower the solar 3 He abundance, that reduction in conjunction with chemical evolution models 
may. We feel justified in claiming that any apparent "crises" in big bang nucleosynthesis is 
rather a (potential) problem for chemical evolution and/or stellar evolution. 
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Figure Captions 



Figure 1: The dependence of the metallicity produced as a function of the deuterium 
destruction factor D/D p , for a large sample of models. The metallicity is 
plotted in solar units. The SFR used was ip oc M gas where the constant of 
proportionality ranges from 0.01 to 1.0 The circles correspond to the choice 
of stellar yields from Maeder (1992) while the crosses correspond to the yields 
of Woosley & Weaver (1993). A power law IMF was chosen with a slope of 
-2.7 for the upper two sets of points and -3.0 for the lower two sets. 

Figure 2: The present day mass function of our adopted model as compared with the 
data from Scalo (1986). 

Figure 3: The evolution of D/H and 3 He/H as a function of time, for a standard model 
of galactic chemical evolution (solid line) and for one which favors massive 
stars early and includes metal enriched outflow (dashed line). Also shown are 
the values of these ratios at the time of formation of the sun, t xs 9.4 Gyr, and 
today, for D/H (open squares) and 3 He/H (filled circles). The present day 3 He 
abundance simply shows the range of observed values; the data point does 
not represent an average. The models were chosen so that D/H is destroyed 
by a total factor of 5, to the present. 

Figure 4: As in Figure 3, for a standard model (solid) and for one in which the stellar 
yields of 3 He at low masses have been reduced (dotted). Deuterium is the 
same in both cases. 

Figure 5: As in Figure 3, for the model with enriched outflow from Figure 3 (dashed) 
and for one with outflow in which the stellar yields of 3 He at low masses have 
been reduced (dotted). 
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